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ABSTRACT 

High resolution blue continuum filtergrams from Hinode are employed to study the umbral fine 
structure of a regular unipolar sunspot. The removal of scattered light from the images increases the 
rms contrast by a factor of 1.45 on average. Improvement in image contrast renders identification of 
short filamentary structures resembling penumbrae that are well separated from the umbra-penumbra 
boundary and comprise bright filaments/grains flanking dark filaments. Such fine structures were 
recently detected from ground based telescopes and have now been observed with Hinode. A multi- 
level tracking algorithm was used to identify umbral dots in both the uncorrected and corrected images 
and to track them in time. The distribution of the values describing the photometric and geometric 
properties of umbral dots are more easily affected by the presence of stray light while it is less severe in 
the case of kinematic properties. Statistically, umbral dots exhibit a peak intensity, effective diameter, 
lifetime, horizontal speed and a trajectory length of 0.29/ QS , 272 km, 8.4 min, 0.45 km s _1 and 221 
km respectively. The 2 hr 20 min time sequence depicts several locations where umbral dots tend 
to appear and disappear repeatedly with various time intervals. The correction for scattered light in 
the Hinode filtergrams facilitates photometry of umbral fine structure which can be related to results 
obtained from larger telescopes and numerical simulations. 

Subject headings: Sun sunspots — umbra — techniques photometric 



1. INTRODUCTION 

The dark umbral background is populated by small, 
bright features called umbral dots (UDs). The 
size of UDs ranges fro m (K'8 to abo ut 0^2 based 
on previous studies by ISobotka et ahl (|l997aj ) and 
Tritschler & Sch midt! (120021) . while a recent work by 
iRiethmiiller et al.l (|2008bf ) shows that the size distri- 
bution of UD diameters is a maximum around 0'/3 or 
225 km, suggest ing that most of the UDs are spa- 
tially resolved. ISobotka et al.l (|l997af ) also observed 
that the larger, long-lived UDs are seen in regions of 
enhanced umbral background intensity. The darkest 
parts of the umbral core, referred to as dark nuclei, 
are often devoid of UDs. Based on their relative loca- 
tion, UDs can be classified as "central" and "periph- 
eral". While the former are seen in the inner regions 
of the umbra, the latter dominate the umbra-penumbra 
boundary. Peripheral UDs are usually brighter than 
the central ones. The intensity of UDs ranges from 
about 0.2 to 0.7 times the normal photospheric in- 
tensity at visible wav elengths. The typical speeds of 
UDs are g 400 m s" 1 (jSobotka et aDl997U iKitai et al. 
20071: IRiethmiiller et al.l l2Q08bt ISobotka fc Puschmann 



m s 1 . UDs do not have a typical lifeti me, with val- 
ues ran ging from 10 min to 2.5 min (jSobotka et al.l 
Il997at IRiethmiiller et "all I2QQ8U lHamedivafal 120111 : 
IWatanabe et al.l l2009h . The spread in the values arises 
from the manner in which UDs can be grouped based on 
t heir size , lifeti me and spatial location. 

IParkerl (|1979[ ) proposed that UDs are manifestations of 
hot non-magnetized plasma pushing its way in the gappy 
umbral field. While the detection of such weak fields in 
the umbra remains elusive, o bservations indicate a reduc- 
tion o f 300-500 G in UDs (ISchmidt fc Balthasar et al.l 
119941 : ISocas-Navarro et al.l 120041 ) with the contrasted 
ones residing in locations where the m agnetic field is ~ 
2000 G and is inclined more than 30° ([Watanabe et al.l 
2009). Central and peripheral UDs exhibit an en- 
hancement i n temperature of 550 K and 570 K 
respectively (Rie thmuller et aTl l2008al ). The measure- 
ments of Doppler velocities in UDs show that peripheral 
UDs h ave an upflow of ~ 0.4-0. 8 km s" 1 (Rimmele 
20041: IRiethmiiller et al.l l2QQ8al ISobotka fe Jurcak 
2009aj ). Central UDs on the other hand exhib it very 



2009). Most mobile UDs emerge near the umbra- 



penumbra boundary and move tow ards the center of 
the umbra ([Riethmiiller et aT1l2008b[ ) with speeds of 700 
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weak downflows (ISo botka fe Puschmann! l20Q9h while 
lHartkorn fc Rimmmeld (120031) dete cted downflows 
of upto 0.3 km s" 1 . lOrtiz et al.l (|201Q[ ) reported 
downflows at the edge of UDs measuring 400 to 1000 
m s" 1 at a s patial resolution o f / /14. 3D MHD sim- 
ulations of Schiissler & Vogler (2006), which model 
UDs as narrow upflowing plumes in regions of in- 
tense magnetic fields, predicted a central dark lane 
in UDs which was subsequently verifie d from high 
resolution ground based ob s ervations (jRimmeld 2008; 
ISobotka fc Puschmann! 120091: lOrtiz et al.l l2010[ Tas well 
as from space ([Bharti et al.l 20071 ). Recent observations 
from the 1.6 m NST indicate that UDs are not perfectly 
circular but possess a mean ec centricity of 0.74 in the 
photosphere (jKilcik et al1l2012l ). 
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The extent to which UDs can be resolved de- 
pends on the spatial res olutio n which is p resently / /14 
(jSobotka fc Puschmannl 120091 : lOrtiz et al.ll2Q10h for a 1 
m ground based telescope at 450 nm. This has been 
made possible with the aid of adaptive optics and post- 
processing techniques to minimize th e contribution from 
"seeing" . On the other hand, Hinode ([Kosugi et al. 2007) 
with a 50 cm aperture in space has a resolution of 0'/23 
at the same wavelength which can operate for long peri- 
ods in the absence of the Earth's atmospheric turbulence. 
Furthermore, since UDs are present in the darkest regions 
of sunspots, contamination by stray light can strongly 
influence photometric investigation of these structures. 
The motivation of this paper is to study the influence of 
stray light on the properties of UDs using Hinode data 
and how its removal compares with existing high resolu- 
tion ground-based observations as well as numerical sim- 
ulations. This exercise has been done taking into account 
the trade off between marginally coarser spatial resolu- 
tion and the uninterrupted "seeing-free" time sequence 
between Hinode and its 1 m ground based counterpart. 
The rest of the paper is organized as follows. The ob- 
servations are described in Section [2] and the results are 
presented in Section [3l In Section |4] we summarize our 
findings and discuss their implications. 

2. OBSERVATIONS AND DATA PROCESSING 

We utilize high resolution blue continuum filtergrams 
of the sunspot in NOAA AR 10944 acquired by the BFI 
(Broadban d Filter Imager) of the SOT (Solar Optical 
Telescope; iTsuneta et al l l2008h , on board Hinode from 
00:14 - 02:34 UT on March 1, 2007. The sunspot was 
located very close to disc center (N0.7W4) at a heliocen- 
tric angle of 6 = 4°. The 1024x512 filtergrams had a 
pixel sampling of / /054 and were taken at a cadence of 
6 s with an exposure time of 102 ms. Initial process- 
ing of Level-0 data to Level- 1 included dark correction, 
flat fielding and removal of bad pixels, and was carried 
out using the "fgjprep" routine in Solar Soft. The images 
were subsequently co-aligned using a 2D cross correlation 
routine. 

In order to reduce the noise in the dark regions of the 
umbra, four successive filtergrams were added to yield a 
sequence of 326 images with a cadence of ~ 26 s. The 
averaging was carried out to kee p the temporal r esolu - 
tion identical to the analysis of IWatanabe et all (2009) 
who had previously studied the same AR with the above 
data set. Applying the running average to the filtergrams 
leads to a reduction of power especially in the high spatial 
frequency domain as shown in the top panel of Figure [T] 
which depicts the azimuthally averaged power spectrum. 
The bottom panel of Figure [1] illustrates the difference 
in power between the single and mean image as a func- 
tion of the spatial scale. Since the power is proportional 
to the square of the intensity, the relative intensity dif- 
ference between the images varies from 2% to 5% for 
spatial scales of 0'/2 and 0'/5, respectively. This implies 
that averaging successive filtergrams should not signifi- 
cantly affect the detection of small scale structures within 
the umbra. We estimate the intensity fluctuations in the 
umbra to be 1.5% of the quiet Sun (QS) intensity. This 
value corresponds to three times the standard deviation 
of a small region in the umbral dark core. Averaging the 
filtergrams reduces the noise by «6%. 
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Spatial Scale [arcsec] 
Fig. 1. — Top: Azimuthally averaged power spectra of single 
(solid line) and average of 4 blue continuum filtergrams (cross sym- 
bol). Bottom: Difference in power spectra of a single and average 
nltergram as a function of spatial scale (grey plus symbols). The 
black cross symbols represent average values of the difference in 
bins of 0"05 with the vertical bars denoting the rras value. 

Insrumental Stray Light: The blue continuum fil- 
tergrams were corrected for instrumental stray light 
using the PSF (Poin t Spread Function) described by 
IMathew et all (|2009f ). This PSF was derived from 
transit observ a tions of Mercury on November 8, 2006. 
IMathew et al.l (j2QQ9h showed that the removal of stray 
light renders an improvement in the contrast of bright 
points in the quiet Sun by a factor of 2.4-2.75 at 430 
nm. As the transit observations were taken close to disc 
center, the same PSF was utilized for removing stray 
light in the blue continuum filtergrams. The PSF is a 
w eighted linear combi nation of 4 Gaussians (see Table 1 
of IMathew et al.l l2009). The deconvolution is carried out 
using an IDL ma ximum likelihood routine 6 ( Rich ardsonl 
1972: lLucvlfT97l . The method uses the instrument PSF 
to iteratively update the current estimate of the image by 
the product of the previous deconvolution and the cor- 
relation between re-convolution of the subsequent image 
and the PSF. The algorithm can be expressed as follows: 

7 = 0*P (1) 

The Image (/) is a result of the convolution between the 
Object (O) and the Point Spread Function P. Given / 

6 Called Max_Likelihood.pro, from the AstroLib package 
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and P, the most likely O can be iteratively determined 
as: 



O t+1 = O* • (J; * conjpj 



(2) 



where C = 0^*P under the assumption of Poisson statis- 
tics. Here t refers to the iteration cycle. 

The corrected and uncorrected sequences were subse- 
quently normalized to the QS intensity. Figure [2] shows 
a typical scatter plot of intensities in a corrected filter- 
gram and the corresponding uncorrected image. Follow- 
ing stray light correction, the minimum umbral inten- 
sity corresponding to the dark umbral core, reduces from 
0.1/qs to 0.05/qs which is above the estimated noise level 
in the the filtergrams. In addition, the fraction of pix- 
els in the umbra having an intensity less than 0.3/ QS in 
both sets of images is « 93%. The removal of stray light 
decreases the mean umbral intensity from 0.202/ QS to 
0.167/q S while the rms contrast increases from 0.059 to 
0.086 which is a factor of 1.45. The middle and bottom 
panels of Figure [3] show the uncorrected and corrected 
image respectively. The umbral region is extracted from 
the sunspot shown as the contoured region (top panel 
of Figure [3]) . The stray light corrected image exhibits 
various fine scale features which are described in the fol- 
lowing section. Hereafter 'uncorrected' and 'corrected' 
will be referred to as UN and CR respectively 7 . 
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Fig. 2. — Change in umbral intensity after removal of stray 
light. Displayed is a scatter plot of uncorrected (I un corr) and cor- 
rected (Icorr) blue continuum intensity in the umbra. Note that the 
axes are scaled differently. The solid and dashed lines refer to the 
intensity histograms of the corrected and uncorrected filtergrams 
respectively with a binsize of 0.01/qs- A magnified version of the 
trailing halves of the histograms is shown in the inset. 



3. RESULTS 

3.1. Umbral Fine Structure 

(a) Short Filaments: Figure [3] depicts three examples 
of short filaments (shown in dashed circles) near the pe- 
riphery of the umbra that closely resemble penumbral 
filaments. These short structures have varying lengths 
with a dark lane and two adjacent brightenings, and a 

7 The movie of the uncorrected and corrected images is available 
at www .prl.res .in/~ eugene/movie- umbral- dots .wmv 
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Fig. 3. — Improvement in image contrast after removal of stray 
light. Top: Sunspot in NOAA AR 10944. The white contour 
refers to the umbral area chosen for analysis. The dotted white box 
representes a selected field-of-view shown in Figure [4] Middle: 
Magnified view of umbra before stray light correction. Bottom: 
Stray light corrected image. The images in the middle and bot- 
tom panels have been scaled identically. The white dash ed circles 
enclose short filaments that are described in Section |3. II 

bead-like brightening at the tip of the filament facing the 
umbra. The width of the filaments shown in the figure 
range from 165 - 200 km where the latter can be con- 
sidered an upper limit for these structures. T hese dark 
filame n ts are similar to the ones rep orted by Ri mmelel 
(|2008f ): ISobotka Puschmannl ()2QQ9h and in numerical 
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simulations ([Rempel et al.|[2QQ9bh . However, to the best 
of our knowledge this is the first time that such fila- 
ments have been seen in Hinode observations. Their 
morphology is different from the tra ditional da r k lan e 
associated with UDs a s obse rved by Rimmele (2008); 
ISobotka fc Puschmannl (|2QQ9ft and lOrtiz et all (|2010h . 
The filament indicated by the largest circle in the bot- 
tom panel of Figure [3] partially obscures a neighboring 
filament to its right. 




Fig. 4. — Temporal evolution of dark filaments and umbral dots. 
The images correspond to the field-of-view indicated by the dotted 
white box in Figure [3] The frames have a cadence of 100s and 
show an UD traveling from the end of one dark filament F2 to the 
flank of another filament Fl. The black plus and cross symbols 
denote the initial and final position of the UD in the time sequence 
respectively. Each major tickmark corresponds to CK'5. 

Figure S] depicts the temporal evolution of a set of dark 
filaments and bright UDs. One can identify at least two 
main filaments labeled Fl and F2. The time separation 
between individual frames is 100 s. The lifetime of Fl 
and F2 is estimated to be 18 and 10 min respectively. 
These values can be regarded as lower limits since both 
of them were present from the start of the sequence. 
The black plus symbol marks the position of a bright 
grain/UD which is seen during the chosen sequence. The 
UD starts at the end of F2 and ends adjacent to Fl on 
its left. The total displacement of the UD from frame 1 
to 11 (cross symbol) is ~0'/15 for the 18 min sequence. 
In general, these localized brightenings/UDs move from 
the tip/end of one filament onto the flank of another. 
This motion is continued till they can no longer be dis- 
tinguished from the background. The dark filaments usu- 
ally maintain their form during the motion of the UDs 
after which they can either diffuse or break up into even 
smaller dark segments. 

(b) Light Bridges: In addition to UDs, sunspot um- 
brae often exhibit light bridges (LBs). These can be 
broadly cl assified into two ca t egorie s, namely - strong 
and faint (jSobotka et al.lll993l fl994h . While the former 
split the umbra into individual cores and represent an 
abrupt change in the umbral morphology, the latter are 
usually less than 1" in width and are composed of a chain 
of UDs. The top panel of Figure [3] shows the umbra to 
be devoid of any large scale structuring, but Figure [5] in- 
dicates a faint LB near the right-hand side of the umbra, 
that is nearly horizontally orientated above an umbral 
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Fig. 5. — Average image constructed from 120 frames in the stray 
light corrected sequence with the white dashed lines indicating a 
faint light bridge. 

dark core. The image is an average of 120 filtergrams 
covering a 50 min duration. The LB consists of several 
bright grains resembling UDs whose width is and 
is close to the resolution limit of Hinode. The average 
grain spacing on the LB is estimated to be 0'/25 while 
the length of the LB is nearly 3.8". 

3.2. Properties of UDs 

This section describes the physical properties of UDs 
which were determined from the time sequence of the 
uncorrected as well as corrected blue continuum image 
sequences. Identifying and tracking of UDs involves the 
following steps- 

Orginal Image Background Level 1 




Level 2 Level 3 Objects Detected 



# 



Fig. 6. — Depiction of MLT algorithm for 3 intensity levels. The 
original image comprises of an arbitrary distribution of intensities. 
The maximum and minimum intensity range is split into 3 inter- 
mediate levels and the objects are sequentially tagged from the 
highest to the lowest level. The object boundary is determined us- 
ing the background image from which the features can be isolated 
as shown in the bottom right panel. 



1. Defining the umbra-penumbra boundary: The se- 
quence of corrected images are added to obtain 
a mean image. After smoothing the mean image 
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using a 11x11 pixel boxcar, the umbra-penumbra 
boundary is defined by a single continuous contour 
corresponding to an intensity of 0.3/ QS . This con- 
tour is used to construct a binary mask which al- 
lowed us to extract the umbral region from individ- 
ual images. 

The identification of UDs was ca rried out using a 
2D MLT (Multi Level Tracking; iBovelet fc Wiehrl 
2001 ) algorithm t hat has been described by 



Riethm uller et al.l (|2QQ8bh for detecting umbral 



dots. The algorithm identifies objects at different 
intensity levels starting from the highest level and 
tagging them uniquely while progressing to lower 
intensities, till the minimum level is reached. The 
number of objects (NOs) detected depends on the 
number of levels (NLs) defined. The latter was 
chosen by implementing the algorithm for different 
levels and counting the total NOs detected. For 
both the UN and CR images, the NOs increase ex- 
ponentially with the NLs. A non-linear least square 
fit provided an optimal value of 33 and 25 intensity 
levels for the UN and CR images respectively. The 
above values correspond to the knee of the best fit. 

3. Defining the backg round image: The thin plate 
spline technique ([Barrodale et al.l 1 19931 ) was em- 
ployed to construct the background umbral image, 
i.e. the intensity distribution in the absence of 
UDs. Each UD is defined by the intensity contour 
corresponding to (i max + / bg )/2, where J max and J bg 
refer to the maximum/peak and background inten- 
sity respectively. Figure [6] depicts the functionality 
of the algorithm based on 3 levels for an arbitrary 
distribution of features. For each UD, the follow- 
ing quantities are determined: 7 max , J bg , / mea n and 
D eii . The effective diameter D eii expressed in km, is 
calculated as ^AA/tt where A is the total number 
of pixels. In addition, the spatial location of the 
maximum intensity (X p , Y p ) is also noted. This 
information is required to track the UDs in the im- 
age sequence. Before saving the above information 
for each UD, the routine also verifies if the UDs 
are separated from the umbra-penumbra boundary. 
Only those features which are at least 2 pixels in- 
ward from the edge of the umbral mask are saved 
and considered for analysis. The fraction of objects 
that did not get filtered using the above criteria is 
less than 5%. 



4. Tracking the UDs in time: Each UD is identified 

in the successive frame if X 1 , Y % and X 2+1 , Y l+1 

p* p p ■ p 

are at the most 1 pixel apart (where i refers to the 
frame index). The condition of 1 pixel separation 
is to ensure that the horizontal speed of the UD 
does not exceed 1.5 km s _1 which is obtained from 
the spatial sampling of 0'/05 and a cadence of 26 s. 
If an UD cannot be identified in the current frame 
the tracking is extended to 2 successive frames. If 
this also fails, the tracking is terminated. 

5. Once tracking is completed, the following addi- 
tional properties are determined from the UD tra- 
jectories: lifetime (T), horizontal speed (V), birth- 
death distance (L hd ) and trajectory length (L ti ). 



The effective diameter as well as the mean, maxi- 
mum and ratio of maximum-to-background inten- 
sities are averaged over the trajectory of the UD. 
The horizontal velocity is calculated as the ratio of 
the trajectory length and the lifetime. 

6. Categorizing UDs into peripheral and central is car- 
ried out using a second boundary that lies «0'/8 
inward from the original umbral mask. UDs which 
originate in between these two contours are labeled 
as peripheral and the rest as central. 




Fig. 7. — Histogram of various UD properties for both time 
sequences. Mean Intensity - /mean, Maximum Intensity - /max, 
maximum-to-background intensity ratio - / m ax//b g , effective Diam- 
eter - D e ff, horizontal speed - V, lifetime - T, birth-death distance 
- L b d and trajectory length - L t j. Bin sizes are 0.05/qs, 0.05/qs, 
0.15, 40 km, 100 m s _1 , 2 min, 50 km, 75 km respectively. The grey 
shaded and unshaded histograms corresponds to the uncorrected 
and corrected time sequences respectively. The y-axis represents 
the fraction of UDs in each bin. 

Figure [71 shows the histogram of the various physical 
properties of all UDs (central and peripheral) whose 
lifetime exceeds 50 s (2 frames). These UDs constitute 
~95% of the total number of objects detected and 
tracked in both time sequences. The histograms have 
been normalized to unity for comparison. The bin size 
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of the quantities corresponding to the UN and CR data 
is identical. We now turn to the quantitative differences 
between the physical properties of UDs derived from the 
uncorrected and corrected sets and whose mean values 
are summarized in Table [TJ 

(a) Intensity: Removal of stray light produces an ex- 
tended tail in the histogram while the peak and minimum 
shift to moderately lower values similar to what was seen 
in Figure El This is more evident in the histogram of peak 
intensity (7 max ). The average value of I mean and 7 max are 
0.24 and 0.26/ QS respectively for the UN set while the 
same for the CR sequence are 0.24 and 0.29/ QS respec- 
tively. The histogram of i mea n peaks at 0.21 and 0.13/ QS 
for the UN and CR data respectively while for J max the 
distribution peaks at 0.16 and 0.19/ QR re spectively which 
is in agreement with lHamedivafal (|201l[ ) . The maximum 
intensity of peripheral UDs is on an average 36% greater 
than that of the central ones after stray light correction. 

The distribution of i max // b g peaks at 1.18 and 1.58 
for the UN and CR sequences respectively. In case of 
the latter the histogram exhibits a conspicuous tail with 
^max being ^4 times greater than J bg for nearly 3% of 
the UD population. As a result, the estimated mean 
value of the ratio become 2.5 in comparison to 1.48 in 
the uncorrected sequence. Such a trend is similar to 
the area histogram of continuum intensit y at 630 nm 
obtai ned for numerically simulated UDs ([Bharti et al.l 
[20101 ). 

(b) Size: Stray light removal reduces the mean effective 
UD diameter from 295 km to 272 km. The histogram 
for the latter is narrower and shows a strong peak at 257 
km while the former is more broader with a maximum 
around 288 km. For the corrected sequence, central and 
peripheral UDs have a mean effective diameter of 298 
km and 239 k m respectively wh i ch is c onsistent with 
the findings of IRiethmuller et al.l (|2QQ8b[ i. The sizes of 
peripheral UDs are smaller than those of cen t ral UD s 
which is sim il ar to that of IRiethmuller et al.l ([20 08b). 
iBharti et al.l (|201Q[ ) state that the average diameter 
of simulated UDs is « 320 km with the histogram 
suggesting that UDs do not have a typical size. In com- 
parison, the histograms of the effecti ve diameter shown 
i n Fi gure and that obtained by Riethmi iller et all 
(2008b) tend to be nearly symmetrical. 
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Relation between L hd and L tj . The scatter plot be- 



tween birth-death distance and trajectory length shown as grey 
cross symbols. The square symbols represent maximum L hd for a 
given L t j within a bin of 100 km. The straight line is a linear fit 
to the scatter of the square symbols. 




arcsec 

Fig. 9. — Trajectories of UDs. The trajectory of UDs is overlaid 
on the time average image from the stray light corrected sequence. 
Blue and yellow represent birth and death respectively. Only UDs 
whose lifetime exceeds 10 min have been shown in the figure. DC1 
and DC2 refer to the two strong dark cores in the umbra. 



(c) Horizontal Speed: Central as well as peripheral UDs 
tend to be mobile with speeds of ^510 and 460 m s _1 
respectively for the uncorrected set. In the corrected 
sequence the speeds of central and peripheral UDs tend 
to be nearly the same, 460 and 440 m s -1 , respec- 
tively, whi ch is in good ag reement with iKilcik et al.l 
(|2012D . IRiethmuller et al.l (|2008b[ ) find peripheral 
UDs to be faste r than central ones by 40 m s _1 while 
IWatanabe et al.l (j2009) report a difference of 170 m s . 
Although the above values show peripheral UDs to be 
slower than central UDs, the maximum speed of the 
former is nearly 200 m s _1 greater than the latter. In 
addition, the difference in the value is small compared 
to the spread in the distribution. The reduced speed of 
peripheral UDs in our case could be attributed to the 
umbral region defined by the mask described earlier. 
Nearly 15% of the population tend to be stationary in 



both sequences. Figure [7] shows that the overall trend of 
the histograms for the UN and CR sets are quite similar 
with some minor differences in the leading half of the 
distribution. 

(d) Lifetimes: UDs do not exhibit a typical lifetime 
as is evident from Figure [7] (second column third 
row). The distribution is exponential with the peak 
coinciding with the smallest lifetime bin. With a cutoff 
of 75 s, our analysis puts the mean value of central 
and peripheral UDs at 8.9 and 7.8 min respectively 
for the corrected image sequ ence. This is cons istent 
with lifetimes ob t ained by lHamedivafal (120081) and 
IRiethmuller et al.l (j2008b). According to lHamedivafal 
(|201l[ ). who obtained a similar exponential distribution, 
central and peripheral UDs have a typical half- life 
of 5 and 3 min respectively while IWatanabe et al.l 
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TABLE 1 

Mean value of physical properties of UDs averaged over their trajectories along with their rms value. The numbers in 

THE PARANTHESES DENOTE THE NUMBER OF UDS IN THAT GROUP. 





All UDs 


Central UDs 


Peripheral UDs 


Parameter 


No Stray Corr. 


With Stray Corr. 


No Stray Corr. 


With Stray Corr. 


No Stray Corr. 


With Stray Corr. 




(1690) 


(1949) 


(953) 


(1096) 


(737) 


(853) 


/mean (J QS ) 


0.24±0.07 


0.24±0.13 


0.21±0.05 


0.21±0.09 


0.27±0.08 


0.28±0.16 


4iax C^Qs) 


0.26±0.09 


0.29±0.17 


0.23±0.06 


0.25±0.12 


0.30±0.10 


0.34±0.21 


■^max/^bg 


1.48±0.34 


2.50±1.26 


1.43±0.28 


2.33±0.91 


1.54±0.39 


2.72±1.57 


D e ff (km) 


295±101 


272±68.0 


337±81.0 


298±54.0 


241±98.0 


239±69.0 


V (km s- 1 ) 


0.49±0.35 


0.45±0.32 


0.51±0.36 


0.46±0.32 


0.46±0.33 


0.44±0.32 


T (min) 


8.80±11.9 


8.40±10.5 


7.10±8.70 


8.90±11.7 


10.9±14.7 


7.80±8.70 


L bd (km) 


94.±120 


112±144 


90.±109 


109±143 


100±133 


115±145 


L t j (km) 


240±314 


221±283 


215±271 


233±309 


271±360 


206±245 



(|2QQ9h report moderately smaller lifetimes of 6.5 and 
7.8 min respectively. These value s however are much 
smaller than the ones obtained by iBharti et al.l (|2010l ) 
who report mean values of ~ 25 min for simulated 
UDs. UDs which are present at the begininng or at 
the end of the time series constitute only 5% of the 
total number and do not influence the lifetime histogram. 

(e) Trajectory: The bottom panels of Figure[7]correspond 
to the birth-death distance and trajectory length of UDs 
respectively. The overall distribution of the above quan- 
tities appear quite similar in both image sequences. 

Nearly 55% of the UD population have trajectory 
lengths of less than 150 km and mean horizontal speeds 
of 375 m s _1 . By comparison, the rest of the group are 
relatively mobile with speeds of 550 m s _1 . The relation 
between L tj and L hd shown in Figure [5J illustrates that 
for a given trajectory length there exists a range of 
birth-death distances. However, the maximum value 
of L hd varies linearly for any L tj or even when the 
latter is averaged within a certain length range. The 
linear relation is shown for the square symbols which 
correspond to the maximum birth-death distance for a 
given trajectory length within bins of 100 km. 



DC1. 



TABLE 2 

Linear correlation coefficients (CCs) between various 
properties for all, central and peripheral uds. the 
values shown in the table correspond to the uncorrected 
sequence while those in parantheses represent 
correlations from the stray light corrected set. ccs 
greater than 0.3 are shown in boldface. 



Correlation 


All UDs 


Central 


Peripheral 


between 








-^max 


- V 


-0.06(-0.004) 


-0.04(0.01) 


-0.01(0.01) 


/max 


- Lhd 


0.24(0.31) 


0.35(0.45) 


0.20(0.25) 


I max 




0.25(0.27) 


0.30(0.39) 


0.22(0.25) 


/max 


- T 


0.33(0.28) 


0.34(0.38) 


0.30(0.29) 


Doff 


- V 


0.21(0.20) 


0.09(0.13) 


0.27(0.26) 


£>eff 


- T 


-0.08(-0.04) 


0.27(-0.13) 


-0.18(-0.06) 


£>eff 


Lbd 


0.21(0.02) 


0.16(-0.16) 


0.31(0.16) 


£>eff 


-L ti 


0.04(0.03) 


0.28(-0.10) 


-0.05(0.08) 


V - 


Lts 


0.25(0.23) 


0.23(0.17) 


0.27(0.33) 


T - 


L^ 


0.86(0.89) 


0.90(0.91) 


0.84(0.86) 


T - 


V 


-0.03(-0.02) 


-0.01(-0.05) 


-0.03(0.02) 



The trajectories of UDs from their point of origin to 
death is shown in Figure [9l These paths have been over- 
laid on the time averaged image from the stray light cor- 
rected sequence. The figure only shows the trajectories 
of those UDs whose lifetime exceeds 10 min. These fea- 
tures constitute 23% of the UD population. Peripheral 
UDs have a tendency to move inward into the umbra with 
a nearly linear trajectory. At least two strong dark cores 
can be identified in the umbra which have been labeled 
DC1 and DC2 in the figure. One finds that UDs tend 
to gather or terminate near the edges of these dark cores 
which is qualitatively in agreement with lWatanabe et al.l 
(2009). Furthermore, there does not appear to be a strict 
segregation of the nature of trajectories on the basis 
of the location/origin of UDs. For instance, UDs with 
nearly linear radial tracks are not necessarily confined 
to the periphery of the umbra while those with squiggly 
trajectories having birth-death distances less than 200 
km are mostly located in the inner regions of the umbra. 
The latter are observed particularly in the bright parts 
of the umbra including the region between DC1 and DC2 
as well as the light bridge on the northern boundary of 



Table [2] lists the linear correlation coefficients (CCs) 
between various properties for both uncorrected and cor- 
rected sequences (values within parantheses for the lat- 
ter). The CCs are calculated between V, T, L hd as well as 
L tj with 7 max and D e ff- CCs greater than 0.3 are shown 
in boldface. The peak intensity of UDs is poorly cor- 
related with its horizontal speed. If one considers the 
correlation between the peak intensity at the point of 
emergence with the average velocity (total distance/total 
time), there is no substantial change from the values cited 
in the table. Brighter, central UDs have larger birth- 
death distances as well as trajectory lengths as indicated 
by the moderate CC which is seen in both the uncor- 
rected and corrected sets. However, the CC worsens in 
the case of peripheral UDs which tends to bring down 
the same for the entire sample of UDs. 

Brighter UDs also tend to be long lived with values 
of 0.38 and 0.29 for central and peripheral UDs, respec- 
tively, as derived from the corrected sequence. A modest 
correlation between size and speed is observed predom- 
inantly in peripheral UDs. In addition larger UDs of 
either kind do not show any preference with lifetimes, 
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birth-death distances and trajectory lengths resulting 
from the corrected images. The same is the case for the 
uncorrected images with the exception of a weak cor- 
relation between size and birth-death distance for pe- 
ripheral UDs. In comparison to central UDs, periph- 
eral UDs exhibit a better, although moderate, correlation 
between mobility and trajectory lengths. The lifetimes 
on the other hand, are strongly correlated with trajec- 
tory lengths. This trend is expected since the trajectory 
length is dependent on the number of frames for which 
the feature is tracked. The lifetime of UDs is poorly 
correlated with its horizontal speed. 

3.3. Spatial and Temporal Coincidence 

The time sequence movie reveals that UDs tend to 
move along trajectories that are traced by subsequent 
UDs. This would suggest that there are regions/pockets 
in the umbra which favor the emergence and disappear- 
ance of UDs. There could also exist a time delay or a 
set of delays between the appearance of one UD and the 
emergence/disappearance of another at the same spatial 
location. This aspect is dealt with in the following man- 
ner. From the information table of each UD, one can ex- 
tract the time interval between the i) birth of two UDs, 
ii) birth of one UD with the death of another and vice- 
versa and iii) death of two UDs, where all three classes 
of events refer to the same spatial location. These events 
can be referred to as birth-birth (BB), birth-death-birth 
(BDB) and death-death (DD) respectively. We confine 
our analysis and discussion to central UDs, since the ori- 
gin of peripheral UDs is predominantly in the pe numbra 
(Sobotka fc Jurca k 2009a; Wa tanabe ena1l2009[ ). 

The top panel of Figure [10] shows the distribution of 
the number of events in the three classes respectively for 
all UDs whose lifetime exceeds 150 s. The histograms 
indicate that there is no unique time lag between succes- 
sive events in either class, although BDB and DD have a 
tendency to have shorter time delays than the BB class. 
Mean time delays of 45.5, 52.0 and 49.8 min were esti- 
mated for the three classes respectively. If the minimum 
lifetime of the UDs is varied in the range 50s - 10 min, the 
mean value of the distribution in all three event classes 
are limited within ±5 min of the above values. Tak- 
ing the average time delay for this range of lifetimes, one 
obtains 46, 53 and 47 min for the each event class respec- 
tively. By comparison, the number of events in BB are 
fewer than the other two for any given minimum lifetime. 
However, for short lived UDs (T <150s), the number of 
BB and DD events are comparable with mean time de- 
lays of 45, 46 and 52 min in the three classes respectively. 

The bottom panel of Figure [T0l shows the spatial distri- 
bution of the occurrences of all three event classes which 
are indicated by the plus, triangle and cross symbols re- 
spectively. The color corresponds to the normalized time 
delay associated with that spatial location for each event 
class. The figure indicates that the pockets of sustained 
emergence as well as disappearance are scattered every- 
where in the umbra with no specific pattern in the dis- 
tribution. The figure also shows that a particular spatial 
location is not limited to the event class as at least 12 
and 20 instances were estimated to be common to BB 
and BDB as well as BDB and DD respectively. While a 
similar characteristic of spatial coincidence of inter event 
classes is seen for UDs with lifetimes of less than 2 min, 



the number of coincidences is higher namely, 19 cases 
each for the above two inter event classes. This would 
suggest that short lived UDs tend to reappear or emerge 
from the point of disappearance of former UDs. 



20 



Birth-Birth 



Birth-Death-Birth 



Death-Death 




20 



40 60 80 100 120 140 
Time delay [min] 




Fig. 10. — Top: Histograms of time delays associated with 
each event class: BB-(grey), BDB-(solid black), DD-(black dot- 
ted). Bottom: Spatial distribution of the occurrences of all three 
event classes (see text) which are indicated by the plus, triangle and 
cross symbols respectively. The color coding indicates the normal- 
ized time delay associated with that spatial location for each event 
class. BB:Birth-Birth, BDB:Birth-Death-Birth, DD:Death-Death. 



4. SUMMARY AND DISCUSSION 

We employ a 2 hr 20 min time sequence of high resolu- 
tion blue continuum filtergrams of the sunspot in NOAA 
AR 10944 from Hinode SOT to determine the properties 
of umbral dots and analyze the umbral fine structure af- 
ter careful removal of instrumental stray light. The mo- 
tivation of our study stems from the fact that stray light 
reduces image contrast which could influence the pho- 
tometric and geometric properties of umbral features. 
After deconvolving the filtergrams with an instrumen- 
tal PSF, whose wings describe the scattered light in the 
telescope, we find that the rms contrast of the images in- 
creases by a factor of 1.45. In addition, the mean umbral 
intensity decreases by « 17%. 
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TABLE 3 

Comparison of a few physical properties obtained previously from Hinode, ground-based observations and numerical 

SIMULATIONS. 



Reference 


Source 


UD Identification Algorithm 


D eff (km) 


-^max /-^bg 


V (km s- 1 ) 


T (min) 


No Stray Light Correction 


Hinode 


MLT 


295±101 


1.48±0.34 


0.49±0.35 


8.8±11.9 


Watanabe et al. (2009) 


Hinode 


Intensity Thresholding 


184 


1.73 


0.44 


7.35 


With Stray Light Correction 


Hinode 


MLT 


272±68 


2.50±1.26 


0.45±0.32 


8.4±10.5 


Hamedivafa (2008) 


0.5 m SVST 


Low-Noise Curvature Detection 


230 




<1.0 


7-10 


Riethmiiller et al. (2008b) 


1.0 m SST 


MLT 


272±53 


1.17±0.1 


0.42±0.20 


10.5±10.5 


Sobotka & Puschmann (2009) 


1.0 m SST 


Low-Noise Curvature Detection 


125 


2.4 


0.34 


9.1 


Kilcik et al. (2012) 


1.6 m NST 


Area & Intensity Thresholding 


254 




0.45 


8.19 


Rempel et al. (2009a) 


Simulations 


Area & Intensity Thresholding 


295 




0.35 


12.9 


Bharti et al. (2010) 


Simulations 


MLT 


420 


2.88 




25.1 



With scattered light accounted for, one is able to iden- 
tify short filamentary structures resembling penumbrae 
but well separated from the umbra-penumbra boundary. 
The features consist of a dark filament with bright fil- 
aments/grains adjacent to it. The width of the dark 
filament is at the resolution limit of Hinode, suggesting 
that these structures are largely unresolved. Their life- 
times are of the order of 10-20 min during which bright 
grains can transit from one filament to another remain- 
ing close to the dark filament. The ex istence of simila r 
features have already been re ported bv iRimme 3(1008); 
ISobotka fc Puschmann! (j2QQ9h based on ground based ob- 
servations from l arger telescopes as w ell as from numer- 
ical simulations (|Rempel et al.ll20"09bft . We however, do 
not find sub- structures in UDs namely, the dark lanes, 
which represent hot convective plumes, since they are 
still below t he detection c apab ility of Hinode. The ob- 
servations of iBharti et al.l (|2007l ) are reminiscent of frag- 
menting light bridges that o f ten exhibit dark cores along 
their axis ([Lites et al.l l2004t iRimmelel 120081 : lLouis et al.l 
I2008h . 

Various physical properties of umbral dots were deter- 
mined after they were identified in individual images and 
tracked in time. To that end, a multi-level tracking algo- 
rithm was applied to both the uncorrected and corrected 
image sequences. The following properties were analyzed 
: mean intensity (J mean ), maximum intensity (I max ), ratio 
of maximum-to-background intensity (I max /I hg ), effective 
diameter (D eif ), horizontal speed (V), lifetime (T), birth- 
death distance (L hd ) and trajectory length (L ti ). The 
mean value of the above quantities derived from the stray 
light corrected image sequence are 0.24/ QS , 0.29/ QS , 2.5, 
272 km, 0.45 km s~\ 8.4 min, 112 km and 221 km re- 
spectively. The distribution of the values describing the 
photometric and geometric properties is more strongly 
affected by the presence of stray light while it is less se- 
vere in the case of the kinematic parameters. 

The time sequence of corrected images was also used 
to determine if umbral dots prefer to appear or disap- 
pear at specific locations in the umbra. Furthermore, 
any specific time delay or a set of time delays associated 
with subsequent emergence/disappearance was also de- 
termined. The events which were sorted as birth-birth, 
birth-death-birth and death-death do not have a unique 
time delay with a majority of events in either group hav- 
ing a value varying between 45-50 min. These values are 
also insensitive to the minimum lifetime of umbral dots. 
The spatial locations where these events occur are scat- 



tered in the umbra and in general do not highlight a spe- 
cific pattern. In addition, there are several pixels where 
events from either class can overlap. This would suggest 
that umbral dots tend to reappear or emerge from the 
location of disappearance of former umbral dots. 

Table [3] gives mean values of a few properties obtained 
previously from Hinode, ground-based observations and 
numerical simulations. The values obtained from the 
uncorrected sequence are in good ag reement with those 
reported by IWatanabe et al.1 (|2Q09j ) for the same data 
set described in this paper, although they employed a 
different identification routine. On the other hand, the 
corresponding values obtained fro m the stray li g ht cor - 
rected images are consistent with [Hamedivafa (2008); 
iRiethmuller et al.l ([2008 oft : iKilcik et all ([2012ft . The ra- 
tio of maximum-to-background intensity (i max // b g) is 
particularly in good ag r eeme nt wit h those obta i ned b y 
ISobotka fc Puschmand (|2QQ9) and IBharti et al.l ([2010ft . 
There appears to be some differences betw een the results 
derived fro m numerical simu lations by Re mpel et al.l 
(j2009aft and IBharti et all ([2010ft which show UDs to be 
larger, slower and long-lived than those retrieved from 
observations. The deviations between the two sets of sim- 
ulations could be attributed to the boundar y conditions 
and treatment of numerical diffusivities ([Kilcik et al.l 
12012ft . In the case of observations, there is general con- 
sensus that UDs are typically 250 km in diameter and 
have lifetimes of 10 min and horizontal speeeds of 450 
m s _1 . Another point of disagreement between simula- 
tions and observations is the ubiquitous presence of UDs 
in the former that do not give any indication of struc- 
turing in the umbra, i.e. faint light bridges and umbral 
dark cores. 

With the exception of intensity vs diameter, none of 
the UD properties are well correlated with each other 
([Kilcik et al.ll2012l ). In simulations, the correlation coef- 
ficient is nearly 0.6 but can vary between 0.02 and 0.3 
depending on the sample of U Ds taken fro m relatively 
bright locations in the umbra ([Kilcik et al.ll2012h . The 
monotonicity seen in simulations, exists with a significant 
scatter between the intensity and diameter of UDs. The 
differences between observations and simulations, needs 
to be addressed with complementary polarimetric obser- 
vations at high spatial resolution (<0 / /15) which could 
provide additional constraints for MHD models. This is 
necessary for determining the physical mechanisms driv- 
ing UDs and t heir kinematics, for wh ich instruments such 
as the CRISP (jScharmer et al.ll2008ft and the GFPI (Gre- 
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gor Fabry-Perot Interferometer; iDenker et al.ll2010l ) will 
be crucial. 
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